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, Abstract 
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In this paper, we show that the exponential integrator scheme both in spatial dis- 
cretization and time discretization for a class of stochastic partial differential equations 
has a unique stationary distribution whenever the stepsize is sufficiently small, and re- 
veal that the weak limit of the law for the exponential integrator scheme is in fact the 
CD ' counterpart for the stochastic partial differential equation considered. 
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1 Introduction 

X 

The convergence and the stability of numerical schemes for finite-dimensional stochastic dif- 
ferential equations (SDEs) have been extensively investigated, see, e.g., Kloeden and Platen 
[T6] and Schurz [20]. Nowadays, numerical approximate schemes for stochastic partial dif- 
ferential equations (SPDEs) are also becoming more and more popular. There is extensive 
literature on strong/weak convergence of approximate solutions for SPDEs. For instance, 
under a dissipative condition, Caraballo and Kloeden [3] showed the pathwise convergence of 
finite-dimensional approximations for a class of reaction-diffusion equations. Applying the 
Malliavin calculus approach, Debussche [6] discussed the error of the Euler scheme applied 
to an SPDE. Greksch and Kloeden [7] investigated the approximation of parabolic SPDEs 
through eigenfunction argument. Gyongy [8], Shardlow [21], and Yoo [23] applied finite dif- 
ferences to approximate the mild solutions of parabolic SPDEs driven by space-time white 
noise. Hausenblas pUl E] utilized spatial discretization and time discretization, including 
implicit Euler, explicit Euler scheme and Crank-Nicholson scheme, to approximate quasi- 
linear evolution equations. Higher order pathwise numerical approximations of SPDEs with 
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additive noise was considered in [13]. For the Taylor approximations of SPDEs, we refer to 
the monograph [T3] . 

However, there are few results on the asymptotic behavior of numerical solutions for 
infinite-dimensional SPDEs although the counterpart for the finite-dimensional case has 
been extensively studied, see, e.g., Schurz [20]. In our present work, we shall investigate the 
asymptotic behavior of certain numerical scheme for a class of SPDEs. To begin with, we 
introduce some notation and thus give the framework of our work. Let (H, (•, •)#, || • \\h) 
be a real separable Hilbert space. Let id# : H — > H be the identity operator, and denote 
(Jif(H), || • ||) and (J£hs{.H), \\ ■ \\hs) by the family of bounded linear operators and Hilbert- 
Schmidt operators from H into H, respectively. In this paper, we consider an SPDE on the 
real separable Hilbert space (H, (•, -) H , \\ ■ \\ H ) in the form 

(1.1) dX(t) = {AX(t) + b(X(t))}dt + a(X(t))dW(t) 

with initial value X(0) = x G H, where W(t) is an //-valued cylindrical id^— Wiener process 
defined on some probability space (f2, J^~,P) with a filtration {^t}t>o satisfying the usual 
conditions, b : H — > H is a Lipschitz continuous mapping, a(x) := cr° + cr 1 (x),x G H, such 
that a G Jf{H) and a 1 : H -> 3f HS (H). 

Throughout the paper we impose the following assumptions: 

(HI) (A,T>(A)) is a self-adjoint operator on H generating an immediately compact Co- 
semigroup {e tA } t >o such that ||e' A || < e~ at for some a > 0. In this case, by [T5l 
Theorem 6.26, p. 185] and [I5j Theorem 6.29, p. 187], —A has discrete spectrum {Aj}j>i 
such that < Ai < A2 < ■ • ■ < Aj < • ■ ■ and lim^oo Aj = 00 with corresponding 
eigenbasis {ej}j>i of H. 

(H2) There exist B x G (0,1) and 8 X G (0, 00) such that JjJ || (-A) ei e sA a°\\ 2 HS ds < 6 X for any 
t > 0, where (—A) 01 := ^2 k>1 A fe 1 (ejj. ® e^) denotes the fractional power of the operator 
-A. 

(H3) There exist Li, L2 > such that 

\\b(x) - b(y)\\ H < L x \\x - y\\ H and ^(x) - cr 1 (y)\\ HS < L 2 \\x - y\\ H , x,y G H. 

(H4) There exists 7 G M such that 

2(x-y, b(x) - b{y)) H + \\a\x) - a\y) f HS < -j\\x -y\\ 2 H , x,y G H. 

By [H Theorem 5.3.1, p. 66], we know that (H1)-(H3) imply the existence and the unique- 
ness of the mild solution to (11.11) . i.e., there exists a unique if- valued adapted process X x (t) 
with the initial value x G H such that 

(1.2) X x (t)=e tA x+ f e^- s)A b(X x (s))ds + f e^ A a(X x (s))dW(s). 

Jo Jo 
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Remark 1.1. In fact, under (HI), (H3) and J * ||e ,sA cr ||^ 5 ds < 62 for any t > and some 
5 2 > 0, (11.11) also admits a unique mild solution on H. While (H2) is just imposed for 
the later numerical analysis. Let a = id Hl and Ax := <9fx for x G V(A) := # 2 (0,tt) n 
Hq(0,7t). Then A is a self-adjoint negative operator and Ae^ = —k 2 ek, k G N, where 
e k(0 '■— (2/71-) 1 / 2 sin £ G [0,7r], A; G N. A simple computation shows that 



^ k =i Jo Z k=i 

Then (H2) holds with 6t = \ E£L# 2 ) 2ei_1 for #1 e (0, 1/4). 
Remark 1.2. By (H3), it is readily to see that 

(1.3) llfe^ll^ + ll^^llL^^l + llxl^), xeH, 

where L := 2((L 2 + L 2 ) V /i) with // := ||6(0)||^ + H^O)!! 2 ^. Moreover, by (H4) one has 

2(x, b{x)) H + \\a\x)f HS = 2(x, b(x) - b(0)) H + \\a\x) - ^(O)^ 



(1-4) + 2(x, b(0)) H + 2(a\x) - ^(0), ^(0))^ + W^Whs 

<-(7-e)||x|| 2 , + 2(L 2 + l + e) yU e- 1 , £6(0,1), x G H, 

where (T,S) H s ■= Y^=ii Te h Se i)H for S,T G 3? H s(H). 

Before establishing the numerical scheme, we further need to introduce some notation. 
For any n G N, let 7r n : H — > H n := span{ei,--- ,e„} be the orthogonal projection, i.e., 
n».x = Y^=i( x ^ e i)Hei,x G H, A n := n n A G J£(H n ),b n := ix n b : H n ->> iJ n and cr n := 7r n cr : 
H n — >■ ^Hs{.H n ). Moreover, throughout the paper, let x n := 7r n x for arbitrary 2 G U, where 
U is a bounded subset of H. 

Consider finite-dimensional approximation associated with (11. ip on H n ~ W 1 



dX n {t) = {A n X n (t) + b n {X n {t))}dt + a n {X n {t))dW{t), 
X n (0) = x n . 



(1.5) 

The spatial approximation (11.51) is also called the Galerkin approximation of (II. ip . Due to 

n n 

n n Ax = ■k ti a( y^(g, ei) H ei) = -^(x, ej) H Xei, x G H n , 



i=l i=l 



it follows that 

(1.6) A n x = Ax, e tAn x = e tA x and (x, b n {y)) H = {x, b{y)) H 

for all x, y G i? n . By (H3) and the property of the projection operator ir n , we have 



\\A n (x -y) + b n (x) - b n {y)f H + \\a\{x) - a\{y)f HS 

< 2\\A n (x - y)\\ 2 H + 2||6 n (x) - K{v)f H + \Wfa) 

< 2\ 2 Jx -y\\ 2 H + 2\\b(x) - b(y)\\ 2 H + \\a\x) - a" 

< 2(A 2 +L 2 1 + L 2 2 )\\x-y\\ 2 H , x,y G H n . 



HS 

HS 



Hence, under (HI) and (H3), (jl.5p admits a unique strong solution {X Xn (t)} t > with the 
starting point x n £ H n . 

Next we introduce a time-discretization scheme for (11.51) . For a stepsize A £ (0, 1) and 

— n A 

each integer fc > 0, compute the discrete Exponential Integrator (EI) scheme Y x ' (kA) « 
X" n (A;A) by setting Y^' (0) := a; n and forming 

(1.7) Y n x f((k + 1)A) := e^iYlfikA) + k(P£ A )A + * n (?^(kA))AW k h 

where AW^ ■— W((k + 1) A) — W(kA), and define the continuous EI scheme associated with 
dL5])by 



;i.8) 



Y£ A (t): = e tA "x n + / e^ A -b n (Y^([s\))ds 
Jo 

+ f e^ A -a n (Y x ^ A ([s\))dW(s) 
Jo 

= e tA x n + f ' e^ A b n (Y:f([s\))ds 
Jo 

+ f e {t - [si)A *n(Y: : A (ls\))dW(s) 



due to p.6p . where \t\ := [t/A]A with [t/A] standing for the integer part of t/A. It is easy 
to see from (11.81) that 

Y x n , A (t) = e^ A Y x ^ A (s) + f e^ A b n {Y^{[r\))dr 

(1-9) ft Js 

+ / e^-^ A a n (Y^ A ([r\))dW(r), 0<s<t. 

J s 

By Y£ A (0) = F^ A (0), we deduce from ([LTD an d (HD that Y" n ' A (£;A) = Y^f(kA), i.e., 
y^* ,A (t) coincides with the discrete EI approximate solution at the gridpoints. 

Remark 1.3. For the finite-dimensional SDEs, the discrete Euler-Maruyama (EM) scheme 
and the continuous EM scheme are standard, e.g., [HI p. 113]. While the roots of construct- 
ing the schemes (jl.8p and (11.91) go back to, e.g., [SJ ITT] . 

For the discrete EI scheme (ll.7p . in this paper we are concerned with the following two 
questions: 

• Given n £ N, for what choices of the stepsize A £ (0, 1) does the EI scheme have a 
unique stationary distribution; 

• Will the stationary distribution of the EI scheme converge weakly to some probability 
measure? If so, what's the weak limit probability measure? 
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In what follows, we shall give the positive answers to these two questions one-by-one. 

It is also worth pointing out that, for the finite-dimensional case, Yuan and Mao [23] 
studied the invariant measure of EM numerical solutions for a class of SDEs, and Yevik and 
Zhao [22] discussed by the global attractor approach the existence of stationary distribution 
of EM scheme for SDEs which generate random dynamical systems. Comparing the EI 
scheme (I1.7P with the EM scheme for the finite-dimensional case, e.g., [TH p. 113], we note 
that the explicit EI schemes (11.7j) is based not only on the spatial discretization but also on 
the time discretization. Moreover, in fll.ip . the linear operator A is generally unbounded, 
and the diffusion coefficient is not Hilbert-Schmidt, which leads to be unavailable of the Ito 
formula. Therefore, our approaches are different from those of [221 E3J- What's more, Brehier 
[2] investigated the existence of invariant measure for semi- implicit Euler scheme (in time), 
and discussed the numerical approximation of the invariant measure for a class of parabolic 
SPDEs driven by additive noise, where the drift coefficient is assumed to be bounded. 

The organization of this paper goes as follows: In Section 2, for a give n G N and a suffi- 
ciently small stepsize A G (0, 1), we show that the EI approximate solution {Y^.' n (kA)}k>o,x n eH, 
admits a unique stationary distribution under the properties (PI) and (P2); Section 3 
is devoting to providing some sufficient conditions such that (PI) and (P2) hold; In the 
last section, we reveal that the weak limit of the law for the EI approximate solution 

{Yx'n (^^)}fc>o,x n e-fln i s i n f ac t the counterpart for fll.ip . 

2 Stationary Distribution for the EI Scheme 

For fixed integer n G N, arbitrary integer k > and T G 3§(H n ), define the fc-step transition 
probability kernel for the discrete EI approximate solution Y^.' (kA) by 



Following the argument of that of [241 Theorem 1.2], we deduce that 
n /\ 

Lemma 2.1. {Y x ' (kA)} k > is a homogeneous Markov process. 

We still need to introduce some additional notation and notions. For a real separable 
Hilbert space (K, \\ ■ \\k), let V(K) stand for the collection of all probability measures on K. 
For Pi, P 2 G V(K), define the metric d]L as follows: 



where L := {/ : K R : \f(u) - f{v)\ < \\u - v\\ K and |/(-)| < 1}. 

Remark 2.1. It is known that the weak convergence of probability measures is a metric 
concept, see, e.g., [12], Proposition 2.5, p. 6]. In other words, a sequence of probability 
measures {Pk}k>i G V(K) converges weakly to a probability measure Pq G V{K) if and only 



p£' A (* n ,r) :=P(C(fcA)er). 



(2.1) 




if lim d L (P fc , P ) = 0. 
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— n A 

Definition 2.1. For a given n G N and a given stepsize A, {Y x ' (kA)}k>o, Xn ^H n is said 



to have a stationary distribution 7r n,A G V(H n ) if lim dL(P^' A (x n , -),7r n ' A (-)) = for every 

k— )-oo 



^ A 

Definition 2.2. For a given n G N and a given stepsize A, {YJ (kA)} k > 0jXnGHn is said to 
have Property (PI) if 

sup sup E||F^ A (A;A)||!<oo 
k>o x„eu 

while it is said to have Property (P2) if 

lim sup E\\T^(kA)-Y^(kA)\\% = 0, 

k ^ 00 x n ,y n eU 

where U is a bounded subset of H n . 

Our main result in this section is stated as follows. 

Theorem 2.2. Assume that (PI) and (P2) hold. Then, for a given n G N and a given 
stepsize A, {Y^,' n (kA)}k>o,x n eH n has a unique stationary distribution 7r n,A G V(H n ). 

Proof. For fixed n G N, we note that H n ~ R n is finite-dimensional, and choose a bounded 
subset U C iJ n such that x n , y n G [/. Following the argument to derive [211 Lemma 2.4 and 
Lemma 2.6], we deduce that 

(2.2) lim sup d L (P£' A 0r n ,-),Pfc' A (2/n,-)) = O, 

and that, together with Lemma [2.11 there exists ir n,A G V(H n ) such that 

(2.3) limd L (Pf(0,),^(-)) = 0. 

K— J-OO 

Then the desired assertion follows from (12. 2p . (12. 3 p and the triangle inequality 

B n,A/ \ _n,A/ \\ ^ j /t^A/ \ Trjr^A/n n\ . i /Tn)n,A/ n n _n,A/ 



d L (^ (x n , ■), tt^(-)) < d L (pr -), pr (o, ■)) + wco, o, ^ A (-))- 



□ 



3 Sufficient Conditions for Properties (PI) and (P2) 

To make Theorem 12.21 more applicable, in this section we intend to give some sufficient 
conditions such that (PI) and (P2) hold. In what follows, C > is a generic constant whose 
values may change from line to line. For notational simplicity, let 

Z n ' A {t):= f e {t -^ )A a° n &W(s) and Y^it) := Y^ A {t) - Z n ' A {t). 
Jo 
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Lemma 3.1. Under (H1)-(H3), 

(3.1) nt n : A (t) - Y^ A ([t\)\\ 2 H < ftA(i + nK: A ([t\)\\ 2 H ), t > o, 

where ft := 3{(A 2 +2L) V (2L(1 + || {-A)~^ || 2 5i))}. 
Proof. Observe from (jl.8p that 

(3.2) Y£ A (t) = e tA x n + f e {t - [si)A b n (Y:f([s\))ds + f e^ A al(Y x ^([sj))dW(s). 

Jo Jo 

This further gives 

t n : A (t) = e (t ~ W) ^i A (W) + T e (4 - W)A ^K A (L,J))d S 



+ 



e^ A al(Y x ^ A ([s\))dW(s) 



L'J 



Then, by the Holder inequality, the Ito isometry and (HI), one has 



(3.3) 



nK A (t)-Y^ A ([t\)\\H 

<3{E||(e^W)^-idHfe A (W)ll^ + E f \\b{Y^{\_s\))f H ds 
1 J It} 



+ E 



W\Y^ A 



HS 



ds 



=: 3{h(t) + I 2 (t) + I 3 (t)}. 
Recalling the fundamental inequality 1 — e~ y < y,y > 0, we obtain from (HI) that 



(e^ A -id H )uf H 



(3.4) 



i=i 



H 



-\ n {t-\t\)\2\ 



U 



H 



< (1 -e 

< XlA 2 \\u\\ 2 H , ueH n . 



Thus we arrive at 

(3-5) h(t) < XlA 2 E\\Y: : A ([t\)\\H- 

Note from the Ito isometry, (HI) and (H2) that 



(3.6) 



n\Z nA (t)\\ 2 H = f ||e^W)A e (,-^ a o || ^ ds < r U -A)^^Af^ A a°Jl s ds 
Jo Jo 



<\\(-A)-°i 



01 l|2 



Thus, by (jl.3p and ( 13. 6 p it follows that 

ut) + h{t) < m\\b(Y: : A {[t\))\\ 2 H + w^ovfimwis} 

(3.7) < 2lA{l + E||^ A (LtJ)||^ + E|^ A (LtJ)||^} 

< 2IA{1 + \\(-A)-^\\% + E\\Y: : A (lt\)\\l}. 

As a result, (13. ip follows by substituting (13. 5p and ( 13. 7p into (13. 3p . □ 

Theorem 3.2. Let (H1)-(H4) hold and assume further that 2a+7 > 0. If A < min{l, (2a + 
7 )7(4p?)}, then 

(3.8) supsu P E||y^ A (t)||^<oo, 

t>0 z„£(7 

where p 1 := 2 + (\Ua - 7| 2 /64 + 2L+ \Ua - 7|/8)/?i + 2(1 + f3 1 + A^L) and C/ is a bounded 
subset of H n . Hence Property (PI) holds whenever the stepsize A is sufficiently small. 

Proof. Note that ( \3.2h can be rewritten in the differential form 

(3.9) dY£ A (t) = {AY^it) + e^ A b n (Y x ^ A ([t\))}dt + e^ A a l n (Y x n : A ([t\))dW(t) 
with Kj™' A (0) = x n . For any v > 0, by the Ito formula we derive from ( 13. 9 p and (HI) that 

E(e^||^ A (t)||^) 

< \\x\\ 3 H + E f e-MY^)^ + 2(Y^%),AY X ^%)) H 
Jo 

(3.10) + 2<^ A (s),e( s -W)^ 

<\\x\\' H + E f^ s {-{2 a -u)\\Y: : %s)\\l 
Jo 

+ 2<Y£ A ( S ),e<--L'^ 

Since 

(3 n) \\Y:: A (t)\\ 2 H = iift A (W)iiH + 2(ft A (W),^; A w -Y: : A ([t\)) H 

+ \\y x n : A (t)-Y:f(it\)\\ 2 H1 

and 

K A W,e( 4 -W)^„(^ A (LtJ)))H^ 

= (Y:: A (it\),b(Y:; A ([t\)))H + (Y x n ^(t) -Y: : A ([t\),b(Y x n ; A ([t\)))H 

- (Z n ' A ([t\)MY x n ; A ([t\)))H + (Y£ A (t),(e^ A -id H )b n (Y: : A ([t\))) Hl 



8 



it follows from ( 13 . 1 j) that 

E(e^||^ A (t)||i)< ||x||^ + E ^^{-(2a-^)||^ A (Lsj)||^+||a 1 K A (Lsj))||L 

Jo 



+ 2(Y:f([s\),b(Y x n : A ([s\))) 



11 

- 2(2a - iy)(Y:f([s\),Y:f(s) - Y^([a\)) E 
-(2a-u)\\Y^(s)-Y: : A ([s\)\\ 2 H 

+ 2(Y x ^ A (s), (e^ A - id H )b n (Y x n f([s\))} H }ds. 
This, together with (11.4)1 . yields that 

ECe^ll^WHlr) < |N|| - (2a + 7 - e - u)E f e-|fe A ( \_s\)f H ds 

t 



+ 2{Y: : %s)-Y:f{[s\)xy x n f{[s\))) 1 

-2<^ a (W)),60^ a (W))>j 



o 



+ E / e"{-2(2a - i/)« A (W),^W - >^ a (W)>h 



o 



n,A/„\ i/n,A/l „ IM|2 



(2a-u)\\Y x n f(s)-Y a 



(3.12) 



+ 2(^ A ( S )-^ : A (LsJ),6K A (bJ))Md, 

+ 2E /V(F«> A ( S ), (e'-W* - id H )b n (Y^ A ([s\))) H ds 
Jo 

+ E f e-{2(^ + 1 + e- 1 )^ - 2{Z"> A ( \_s\\ b(Y^ A ( [s\))) H 
Jo 

- 2( 7 - e)(Z^ A ([s\),Y: : A ([s\)) H - ( 7 - e)||Z"' A (L,J)|||}d S 



=: Ji(t) + J 2 (t) + J 3 (*) + J 4 (*). 
By the elemental inequality: 2a6 < ko 2 + b 2 /n, a,&6l,K>0, and ( 13.11) . we arrive at 

•M*) < E f ^{^\\Y:; a {[s\)\\ 2 h + 2- l r- 1 A^||&(K c -^( [aj))lllr 



+ {(|2a - H 2 + 2l)A-s + |2a - H}lfc A (*) - *£ A (kl)lllr}d* 

<E^-{2A^||^ A (L,J)||| + 2- 1 AUAi||Z^(L,J)||| 

+ {(|2a - H 2 + 2l)A-i + \2a - v\}\\Y^{s) - Y^{ \a\)f H )&s, 
where in the last step we have used ( 11.31) . Combining ( 13. ip with ( 13. 6p . we thus obtain that 

Ut) < / <e ^{{2 + (|2a - H 2 + 21+ |2a - H)MA^ft A ( L*J)llfr 



(3.13) 



{1 + \\{-A)- e ^f5 l + (|2a - u\ 2 + 2L+ \2a - u\)^}A^ds. 
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On the other hand, we deduce from (11.31) . (13.41) . (13.61) and (13.111) that 

Mt) <E ^ e^{Z\^ Ltj ) HI, + (?-^( W ), n^(^) - L*J )>^ 

Jo 

+ ^\K A (t)-Y^(ml + ^\\(e^ 

(3.14) Jo 

+ A"* || (e<- ^ A - id H )b n (Y x ^( [s\ )) \\%}ds 

< [ t e-{2(l + (3 l + XlL)A^\\Y x ^([t\)\\ 2 H 
Jo 

+ 2((3 l + XlL(l + \\(-Ar^\\%))A^}ds. 
Furthermore, due to (II. 3p and (13. 6p . for arbitrary k > one has 

Ut) < E /V{2(L* + 1 + e" V" 1 + 2\\Z^{[s\)\\ H \\b{Y:;\ [s\))\\ H 
Jo 

+ 2\j-e\. ||^' A (W)lkll>i' A (W)lk + 17 " e| • ll^ A (W)lllr}d« 

< E fe vs {2{L\ + 1 + e" 1 )^" 1 + ^\\Z^{\_s\)f H + K\\b{Y^{\_s\))\\l 

Jo 

+ | T - e| 2 ^ 1 ||z"" A (L s J)||H + *|fe A (bJ)lli + It- *\ ■ ll^ A (W)l&}d* 

< / e vs {KL + {K- 1 + 2KL+\ 1 -e\ 2 K~ l + \ 1 -e\)\\{-AY^\\ 2 5 l 
Jo 

+ 2{L\ + 1 + e" 1 )^ 1 + (1 + 2L)KE\\Y^ A ([s\)\\ 2 H }ds. 
In particular, taking e = u = (2a + r y)/8 and k = {2a + 7)/(4(l + 2L)) yields that 

(3.15) J 4 (t) < [\^{A-\2a + 1 )E\\Y: : A ([s\)f H + C}ds. 

Jo 

Putting f[3TT31 - f[3TT5D into fl3TT2j) . we deduce that 

n^Kf{t)\\%)< \\x\\h + C feTds 

(3.16) 



2a + 7-2 Pl l 



■E / enft A (W)ll|d S . 

JO 

For A < (2a + 7) 2 /(4p 2 ), it is trivial to see that 2a + 7 - 2piA^ > 0. Thus we have 



sup sup E(||l£*(t)||j,)<oo. 

t>o x„eu 

Finally, flE2i follows by recalling Y^ A (t) = Y^ A (t) - Z nA (t) and (JUS]). □ 
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Theorem 3.3. Let the assumptions of Lemma [3T21 hold. If A < min{l, (2a + l) 2 / (^p^)} , 
then 

(3.17) lim snp E||l£ A (t) - Y£ A (t)\\ 2 H = 0, 

where p 2 := 6(A^ + L)(\2a — j\ + 1) + 3 + 7L + X^L + 6A^ and [/ is a bounded subset of H n . 
Hence Property (P2) holds whenever the stepsize A is sufficiently small. 

Proof. Let 

Note from (JUBJ that 

+ T e (t " W)A (6„.K A ( W)) " 6»(i^ A (W)))d* 
./[tj 

+ T e (t - W)A KK A (W)) - ^(^ A (W)))dW( S ). 
■/[*! 

Following the argument of that of (13. ip . we derive that 

(3.18) n\z:> A yn (t) - z^ yn {\t\)f H < ^i+L)mz:tvSm\v 

For v := (2a + j)/2, by the Ito formula it follows from (TO]) , (HI) and (H4) that 

E(e^||^l(t)|||r) < ||x- * + I /E f e-||Z^ n ( S )||l,d S 

Jo 

+ E jT e-{2(^ n ( a ), (a))* 

+ 2(^1 ( W ), &K A ( W)) " &Q£ A ( W))>* 

+ lk 1 K A (W))-^ 1 K A (W))ll^ 

+ 2(^;l(s), (e'-WM - id H )(6„K' A ( W)) - ^« A ( kl)))M<k 

< \\x - y\\ 2 H - (2a + 7 - i/)E /V||3^( [s\)f H ds 

Jo 

+ E /V{-2(2a - ^(^(W),^^) - 



(2a - 



t 



+ 2{^iW -^KW),*Q£*(W)) -6« A (W))Md* 

+ 2E f e-(^( S ), ( e C-WM - idH)(6„(^; A ( )) - Ml£ A (kl)))>* 
./o 

-:Ji(t) + J 2 (t) + J 3 (t). 
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where we have also used the fl3~TTj) with Y£ A (t) replaced by Z^ A yn (t). By (H3) and (l3~T8j) . 
one has 

3.(*)<e /* £ e-{z\^||^^(L s J)||^ + W)) - &(>X'^(L^J))II^ 

Jo 

+ {\2a-u\ + (\2a -v\ + 1)A^}||Z^( S ) - Z:±([ S \)\\ 2 H }ds 

< {6(A 2 „ + L)(\2a -v\ + l) + 1 + L}A^E f ^\\Z^ yn {[s\)\\ 2 H ds. 

Jo 

On the other hand, carrying out a similar argument to that of (13.141) leads to 

J 3 (t) < 2E f e vs {^\\Z:±{s) - Z^{\s\)f E + A^ A » - Z^([s\), Z^ A n ([s\)) H 
Jo 

< + + + 6L)AiE f e vs \\Z:±{[s\)\\ H &s. 

Jo 

Hence we arrive at 

E(e-||^ n (t)|||) < \\x - yf H - 2a + 7 ~ 2P2A " ^[ e^ll^:l(bJ)H^ 
and then the desired assertion (13.171) follows by A < min{l, (2a + / ~f) 2 / \^p%)} ■ □ 

4 Weak Limit Distribution 

In the previous section, we give some sufficient conditions such that (jl.7p has a unique 
stationary distribution 7r n,A G V(H n ) for a fixed n and a sufficiently small stepsize A G (0, 1). 
In this section we proceed to discuss the weak limit behavior of 7r n,A G V(H n ) and give 
positive answers to the following questions: 

• Will the stationary distribution 7r n,A (-) converge weakly to some probability measure 
in V(H) whenever n — > oo and A — > ? 

• If yes, what is the weak limit probability measure ? 

Denote {X x (t)}t>o jX eH by the mild solution of (11.11) starting from the point x at time 
t = 0, which is a homogenous Markov process. For any subset T C 38(H) and arbitrary 
t > 0, let ¥ t (x,T) := P(X x (t) G T). 

Definition 4.1. {X x (t)} t >o. x eH is s& id to have a stationary distribution ir(-) G V(H) if 
limd L (P t (x,-),vr(-)) = 0. " 

t— >oo 

To reveal the limit behavior of 7r"' A (-), we first give several auxiliary lemmas. 
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Lemma 4.1. Let (H1)-(H4) hold and assume further that 2a + 7 > 0. Then the mild 
solution {X x {t)} t > 0x( z H of (11. ip has a unique stationary distribution ir(-) G V(H). 

Proof. We remark that [U Theorem 3.1] investigates the stationary distribution of (11.11) with 
cr° = 0, i.e., the diffusion coefficient there is a Hilbert-Schmidt operator. For cr° 7^ 0, note 
that a is not Hilbert-Schmidt. Therefore [H Theorem 3.1] is unavailable for (II. ip . Let 

(4.1) Z(t) := [ e {t - s)A a°dW(s) and X x {t) := X x (t) -Z(t). 

Jo 

Then (II. ip can be rewritten in the form 

(4.2) dX x (t) = {AX x (t) + b(X x (t))}dt + a\X x {t))dW{t). 

To be precise, (14.21) is first meant in the mild sense. But under (H1)-(H3) it also has a unique 
variation solution, and therefore the Ito formula applies to H-X^i) ||#. Carrying out similar 
arguments to those of Theorem 13.21 and Theorem 13.31 respectively, for some bounded subset 
U C H we deduce that 

(4.3) supsupE||X x (t)||^ < 00 

t>0 x&U 



and 



lim sup E||X x (t) - X y {t)f H = 0. 



x,yeU 

Then, following the argument of that of [U Theorem 3.1] yields the desired assertion. □ 

Lemma 4.2. Let (HI) and (H2) hold and assume further that there exists 62 > and 
6 2 e (0, 1) such that 

(4.4) / ||e s V||^ 5 ds < 5 2 A e2 . 

Jo 

Then 

(4.5) supE||Z(t) -Z([t\)\\ 2 H < CA eiAe ' 2 , 

t>o 

where C > is a constant independent of A. 

Proof. Recall from [T9l Theorem 6.13, p. 74] that there exists C\ > such that 

(4.6) ||(-A) Ql e^|| < dr a \ \\{-A)- a2 {l - e tA )\\ < C x t a \ 
for arbitrary ol\ > 0, «2 G [0, 1], and that 

(4.7) {-A) aM x= (-A) a *(-A) a4 x, ieD((-A) 7 ), 
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for any 03,04 £ E, where 7 := max{a 3 , 04, 0:3 + 04}. In the light of the independent 
increment of Wiener process and the Ito's isometry, 

n\Z(t)-Z{[t\)\\ 2 H = / LtJ ||(e ( '- LiJ)A -idH)e(W-)V||^d S + f He^VU^d*. 

Jo J\t\ 



This, combining (H2), ( jOJ) , ( Oj) with (TCTD, y ields that 



E||Z(t) - Z(LtJ)||| < / W \\(-A)-^(e^ A - id H )f ■ \\(-A) 6 ^-^ A a \\ 2 HS ds 



A 



\e sA a°\\ 2 HS^ 
pit] 

<C 2 A 261 / \\(-A) ei e sA a°\\ 2 HS ds + 25 2 A 62 
Jo 

< (C 2 1 5 1 + 5 2 )A e ^ e \ 

and therefore the desired assertion follows. □ 

Remark 4.1. Let a = id# and A be the Laplace operator defined in Remark 1 1.1[ A straight- 
forward computation shows that 

/A I 00 -I 

ii^iiL^ = ^EF (1 - e " 2fc2A) - 
fc=l 

Recall that for arbitrary 5 £ (0, 1) and 2, y > 

(4.9) |e- x -e-^| < Ix-yl 5 . 

It then follows from and (jOJ) that 

fc2(l-<5) • 



/A 00 -, 

||e-||^d S <2^A^^i 



Hence, flO} holds with 5 2 = 2 5 - 1 J2T=i ^h) and ^ = *e (0, 1/2). 
Lemma 4.3. Let the assumptions of Lemma [4.11 hold and 

(4.10) r := oT x L x + (2a)~ 1/2 L 2 £ (0, 1). 
Then 

(4.11) supE||X,(f) - Y^{t)\\ 2 H < C{\J^ + A 9 ^}, 
t>o 

where C > is a constant dependent on x £ if but independent of n and A. 
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Proof. By ([O]) and it follows that 

(4.12) supE\\b(X x (t))\\l + supE\\a\X x (t))\\ 2 HS < C. 



t>0 



Note that (E|| • H^) 1 ^ 2 is a norm and recall from 0, Theorem 202] the Minkowski integral 
inequality: 

2n 1/2 /"* / n 1/2 



E 



F(s)ds 



< 



E|F(s) 



ds, t > 0, 



where F : [0, oo) x Q — > R is measurable and locally integrable. Then, applying the Ito 
isometry and using (HI), we obtain from (jl.2p that 

(Eiix^-x.atj)^) 1 / 2 

< \\ e lt]A {e (t-lt»A_ l}x y 



LfJ 



(4.13) 



+ 



( E || e ( L*J s)A | e (t- L*J M _ id H }b{X x (s) ) || ^) 1/2 ds 
W E\\e^^ A {e^~^ A - id H }a 1 (X x (s))\\ 2 HS ds'" 



+ / (E\\b(X x (s))\\ 2 H y/ 2 ds 

J\t\ 



E\\a\X x {s))f HS ds 



1/2 



Let p 
has 



=: Ft{t) + F 2 (t) + F 3 (t) + F 4 (t) + F 5 {t). 
h A 2 )/2. In view of (Oil . (HTTjl . (HI) and the boundedness of (-A) _(1_p/2) , one 



2 

I 7-i 



= ||(-^)-( 1 ^/ 2 ) e W^(-A)-^ 2 {e(*-W) A -id ff }(-A)x| 

< ||(_^)-(i-P/2) e L*JA||2 . ^_ A y P /2 {e (t-[ti)A _ idH} (_ A)x f H 

< C\\{-A)-^- p /^f ■ \\Axf H A p . 

Also, by (USD and (gZD, we obtain from (Q21I that for e (0, 1) 

5 

•L*J 

(-A)^{e ( '" W)A -l}||ds 



(4.14) 



fc=2 

< CA 1/2 + C 
•L*J 



-v4) p e e(L * J " s)A || ■ ||e (1 " 9)(L '^ s ) A | 



.^)P e e(LtJ-)A||2 . || e (l-fl)([tJ-)A||2 . ^_ A y P{e (t-it})A _ 1} ||2 ds 







1/2 



L'J 



C CA 1/2 + CA P / (0sr'e- Q(1_fl)s ds + CA P 



L*J ~ \l/2 



Observe that 
-L'J 



o 



s -p e -a(l-6)s ds < _ / s -p e -^ ds = ( a ( X _ 0))P-lr(l - p), 
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and similarly 

•L*J 



/ S - 2p e~ 2a{1 - 9)s ds < (2a(l - ^)) 2p " 1 r(l - 2p), 
Jo 

where T(-) is the Gamma function. Hence 



k=2 

This, together with the estimate of Fi(t), gives that 



supE\\X x (t) - X x ([t\)\\ 2 H <CA 

t>o 

Noting that X x (t) = X x {t) — Z(t) and utilizing ( 14. 5p . one has 

(4.15) supE\\X x (t) - X x ([t\)\\ 2 H < CA 9 ^. 

t>o 

Since 



\\(id H - 7r n )(-Ay dl u\\ 2 H = I ^ A^(n,e fc ) fl e fc 

fc=n+l 

<A- 2ei H||, ue#, 
we arrive at 

(4.16) \\(id H -7r n )(-A)^\\ 2 <\; 2e \ 

By virtue of the ltd isometry, (H2), ( KIM . (TO and ( ETTD . it follows that 

E\\Z(t) - Z n ^{t)\\ 2 H < 2 f \\e sA (id H - K n y\\ 2 HS ds 

Jo 

+ 2 f ||(-A)^ 1 (id H -e( s "W)^)(_^)ei e (^M a O||^ ds 
(4 1?) < 2||(idn - n n ){-A)^\\ 2 \\(-A^e sA a°\\ 2 HS ds 



+ CA 2 ^ I \\{-Ay^ A alf HS ds 
Jo 

< C(\\(id H - n n )(-Ay ei \\ 2 + A 2e f \\(-A)^e sA a°\\ 2 HS^ 

Jo 



<C{K 2e ^ + A 2 ^). 
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Following the argument of (I4.13p . we have 

(nx x (t)-Y x ^(t)\\%)^ 

< \\e tA (id H - Tr n )x\\ H 

+ f \\^~ s)A {id H - n n )\\(E\\b(X x (s))\\H) 1/2 ds 
Jo 

"* \ 1/2 



+ 



(jf lle^^-vrOfEHa 1 ^.^))!!^)" 

+ / < ||e ( '" s)A ||(E||&„(X x .( S ))-6 n (X :l .(L S J))||^) 1/2 ds 
>/ o 

1/2 



(4.18) 



|e ( *- )A || 2 E|ki(X iB ( 8 ))- i(X,(L a J))||| r 5<ia 

+ / l|e ( '- s)A ||(E||6„(X x .(L S j))-6 n (F x !:' A (L S J))||i) 1/2 d S 
./o 

+ (£||e^|| 2 E||^(X x .(L,J))-^(F™' A (Lsj))||Ld S 
+ f \\e^ A {id H - e^ A }\\ (E||6(^ A ( bJ )) |||) 1/2 d S 



1/2 



, ft \ 1/2 

+ ( ^ lle^id^-eC'-W^II^II^^CW))!!^) 

i=l 

A straightforward computation shows that 

oo 

||e M (id„ - it n )uf H = e ~ 2Al '(^ ueH. 

i=n+l 

This further gives that 

(4.19) ||e* A (id^-7r n )|| 2 < e ~ 2A "' 

and that 

/ _2L e -2A,i .1/2 

(4-20) G x {t) < ( £ ^xT^Mh) < K'WMh 



i=n+l 1 



by recalling that {Xi}i>i is a nondecreasing sequence. By ( 14. 12ft and ( 14. 19ft . one has 

G 2 (*) + G 3 (t) < C [ ||e^ A (id H - 7r n )||d S + C( / ||e^ A (id^ - 7r n )|| 2 d S 
(4.21) Jo ^ 



< C I'e-^-lds + c(^e" 2A "^d S ) 1/2 < C{\~ 1 + A~ 1/2 ). 
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Taking (HI), (H3) and (14.151) into account gives that 

G 4 (t) + G B (t) < CA&W{ jf ||e^||d S + ( jf ||e^|| 2 d S ) V2 } 
< (7A (eiAe2)/2 . 
Next, note from (HI) and (H3) that 
G 6 (t) + G 7 (t) 

< sup(E||6(X x (L S J))-6K A (W))ll^) 1/2 f lle^Uda 
+ snp (na\X x ([s\))-a\Y^ A (ls\))\\ 2 H ) 1/2 ( f ||e ( ^ )A || 2 d S ) ^ 

0<s<t V Jo 7 

(4.23) < a" 1 snp <p\\b(X x (\s\)) - b(Y£*{ [a\))\\ 2 H )^ 

0<s<t 

+ (2ar 1/2 snp (EUa^X.CLaJ)) - a\Y^{ \_s\))f H )^ 

0<s<t 

<r sup (E||X X ( S )-F« : A ( S ) ii 2 ^ 

0<s<< 

<r sup (E||X.( S )-^ A ( S )|| 2 ,) 1 /2 + r sup (Ep(s)-Z^(s)||^ 

0<s<t 0<s<i 

where r G (0, 1) is defined by (I4.10p . Following the argument of (I4.14p leads to 

(4.24) G 8 (t) + G 9 (t) < CA (eiA92)/2 . 
Substituting ( l4T20]) -( l4T24l) into ( 14181) yields that 

bup(E||X,(*) - ?^; A (t)||i) 1/2 < C(A;^ 2 + A^V 2 ) 

t>0 

due to r G (0, 1). Consequently the desired assertion follows from (14.171) . □ 

Theorem 4.4. Assume that (H1)-(H4), (|4~4~1) and (14.101) hold. Then, there exists a A n such 
that lim^oo A n = 

lim d L (7r"< A "(-),7r(-)) = 0. 

n—>oo 

Proof. Fix x G H and let e > be arbitrary. By Lemma 14.31 there exist a sufficiently large 
neN and a A n sufficiently small such that such that 

d L (P fcSn (x,-),P:' S "(x n ,-))<e/3. 

For the previous n G N, by Theorem I2.2[ there exist a sufficiently small A n and 7\ > such 
that 

d L (P^ A "(x n ,-),vr^(-))<e/3 
18 



whenever kA n > T\. Furthermore, due to Lemma [4.11 there exists T2 > such that 

d L (P t (x, •),?!-(•)) < e , t>T 2 . 

Let T := T\ V T 2 , A n = A n A A n and /c = [T/A n ] + 1. Then the desired assertion follows 
from the triangle inequality 

d L (7T"< A (-), 7T(0) < d L (P fcA (x, •), 7T(0) + 4(PfeA(x, Pfc' A (^n, ")) 

+ d L (P^ A (x n ,-),^ A (-))- 

□ 

Remark 4.2. For the finite-dimensional case, finite-time convergence of numerical scheme is 
enough to discuss the limit of stationary distribution of numerical solution [T8| Theorem 
6.23, p. 266]. While for the infinite-dimensional case, we need the uniform convergence of EI 
scheme (11.71) to reveal the limit behavior of 7r n,A , which is quite different from the finite- 
dimensional cases, and therefore (14.101) is imposed. On the other hand, for the finite-time 
convergence of EM scheme (jl.8p . condition (14.101) can be deleted by checking the argument 
of Lemma 14.31 and combining with the Gronwall inequality. 

Remark 4.3. By following the procedure of this paper, numerical approximation of station- 
ary distribution of SPDEs with jumps can also be discussed, which will be reported in 
forthcoming paper. 
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